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The 8 TeV LHC Higgs search data just released indicates the existence of a scalar res- 
onance with mass ~ 125 GeV. We examine the implications of the data reported by AT- 
LAS, CMS and the Tevatron collaborations on understanding the properties of this scalar 
by performing joint fits on its couplings to other Standard Model (SM) particles. We dis- 
cuss and characterize to what degree this resonance has the properties of the SM Higgs, 
and consider what implications can be extracted for New Physics in a (mostly) model- 
independent fashion. We find that, if the Higgs couplings to fermions and weak vector 
bosons are allowed to differ from their standard values, the SM is ~ 2a from the best 
fit point to the current data. Fitting to a possible invisible decay branching ratio, we find 
BRinv ~ 0.05 ± 0.32 (95% C.L.). We also discuss and develop some ways of using the 
data in order to bound or rule out models which modify significantly the properties of this 
scalar resonance, and apply these techniques to the current global dataset. 



I. INTRODUCTION 

Particle physics entered a new era with the announcement of the discovery of a new boson [[T]| 
based on excess events in several Higgs search channels using 7 + 8 TeV LHC data collected in 
201 1-2012. In light of this discovery, it has become obvious that the question of central importance 
to address now is - what are the properties of the scalar field responsible for the observed excesses? 
The answer to this question determines if this field corresponds to the Standard Model (SM) Higgs, 
with the specific SM mechanism of elegantly breaking electroweak (EW) symmetry, or whether 
a more complicated mechanism is involved in EW symmetry breaking. We study this question 
in detail in this paper, characterizing to what degree a SM Higgs is consistent with the current 
global dataset and presenting several results on the properties of the scalar field. Besides updating 
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and expanding our past results flU, we also present new analyses that emphasize the power of the 
growing dataset to bound and rule out alternative models or to give hints of New Physics (NP). 

It is worth emphasizing that it is very important to specify (and justify) a coherent theoretical 
framework in which to study the emerging evidence for the scalar field. However, without knowing 
the ultraviolet (UV) origin of this field, we do not know what effective field theory (EFT), or 
complete model, should be used to fit the data. We emphasize that at this time, the existing 
experimental evidence is not sufficiently strong to directly assume that the scalar resonance is the 
SM Higgs boson, ascribing any deviations in the measured properties of the scalar field directly 
to the effects of NP interactions expressed through higher dimensional operators. Although this 
is certainly one possible interpretation of the data (and we will examine the implications of Higgs 
data for NP in this framework), we emphasize that, in general, one should not assume what one 
wishes to prove. 

Nevertheless, in formulating a theoretical framework for this study, a wealth of other exper- 
imental results that are also sensitive to the properties of scalar fields at the weak scale can be 
distilled into simple physical guiding principles. These are, namely, approximate Minimal Flavour 
Violation (MFV) |I3]-|7]|; respecting the soft Higgs theorems of Refs. |l8l[9l (i.e. the scalar couples 
to the SM fields in proportion to their masses); and an effective breaking of custodial symmetry, 
SU(2)c, lfT0] - fT2l approximately as in the SM. Directly incorporating these principles in the for- 
mulation of the effective Lagrangian allows us to restrict our attention to a few simple cases. In 
order to establish experimentally the properties of the scalar resonance in a model-independent 
way, one can utilize the effective field theory of the chiral EW Lagrangian coupled to a scalar field 
that was emphasized in Refs. [l2l[l3l to study recent Higgs signal- strength data.^ Depending on the 
assumptions of the UV origin of such a Lagrangian, one is lead to various sets of free parameters 
to fit the data when studying the consistency of the SM Higgs hypothesis with the current dataset. 
We discuss and utilize this framework extensively in this paper to examine the properties of the 
scalar field emerging from the data. 

We also emphasize that with the discovery of a new scalar resonance, one can also use the signal 
strength properties of the scalar field to bound and rule out models that provide too few signal 
events as well as models that provide too many signal events. Further, one can also exclude allowed 
parameter space due to the degree of tension within the dataset that depends on the properties of 



For other model-independent approaches to the determination of the Higgs couplings, see Refs. lfT4ll23l . 
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the scalar field. As the dataset evolves, these techniques become complementary to direct fits 
on the signal strength dataset. These bounds can be quantified by excluding parameter space in 
the allowed couplings of the scalar using a Gaussian probability density function approach. We 
develop and apply such an approach in this paper. 

The outline of this paper is as follows. In Section|ll]we discuss the EFT framework we employ, 
and the implicit assumptions about the UV origin of the scalar field that are adopted when fitting 



the data with various free parameters. In Section III we review and discuss the manner in which we 



treat the scalar signal strength data and electroweak precision data (EWPD), while in Section IV 



we present results, based on our fit method, of the status of the SM Higgs hypothesis. In Section 



IV B we discuss some of the implications of Higgs signal strength parameters for beyond the SM 
(BSM) physics, expressed through model-independent free parameters. In Section |V] we discuss 
novel and complementary methods to identify the allowed parameter space, and in Section VI we 
conclude. 

II. THEORETICAL FRAMEWORK 

An effective chiral EW Lagrangian with a nonlinear realization of the SU(2)lxU(1)y symme- 
try gives a minimal description of the (non-scalar) degrees of freedom of the SM consistent with 
the assumptions of SM-like SU(2)c violation and MFV The Goldstone bosons eaten by the W^, Z 
are denoted by vr" (where a = 1,2,3), and are grouped as 

S(x) = 6^"''""/% (1) 

with aa the Pauli matrices and v = 246 GeV. In this approach, the EW scale v, which sets the mass 
of fermions and gauge bosons is introduced directly into the Lagrangian. The S(x) field transforms 
linearly under SU(2)l x SU(2)r as S(x) L S(x) where L, R indicate the transformation on 
the left and right under SU(2)l and SU(2)r, respectively, while SU(2)c is the diagonal subgroup 
ofSU(2)L X SU(2)r. 

Adding a scalar field h to this theory is trivial. One chooses h to transform as a singlet under 
SU(2)c and a derivative expansion of such a theory is given by 

^eff = l{d,hf-V{h) + jTT{D,^^D^J:) 



h 

2a- + b — + b3 — + 

V V 



h 

V f ^ 



Vij Ur 



h.c. • • • , (2) 
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with 

Although we use the notation h, we do not assume in principle that this scalar field is the Higgs, 
or that the scale v is somehow associated with the vacuum expectation value of this field - as this 
is what we seek to establish from the data. As is well known, the a and Cj parameters control the 
couplings of h to gauge bosons and fermions, respectively, and therefore, play a crucial role in 
the phenomenology of single h production. Note that in previous fits, and in the majority of this 
work, the assumption Cj y^-'^ = cy^f" is used and no distinction is made between the rescaling of 
the /i-coupling to the and rf* quarks. We will relax this assumption later on. This Lagrangian 
is common in the study of composite models and its relevance has been emphasized recently in 
Refs. [|24] - l26ll . It is also appropriate to study a pseudo-Goldstone boson (PGB) emerging out of 
an approximately conformal sector [|27ll28ll . or as the low-energy EFT arising in many scenarios 
where the Higgs is a composite PGB that emerges from the breaking of a larger chiral symmetry 
group ||29] - [34ll . We emphasize that this EFT setup can be matched to many UV frameworks and, 
being quite general, we do not confine ourselves to any particular UV scenario.^ 

We approach the data in this way to be as model-independent as possible. However, even 
specifying the free parameters that one will use to fit the data introduces implicit model de- 
pendence. One can nevertheless broadly characterize certain parameter choices. As this is an 
EFT, Ceff is non-renormalizable. Here (and throughout this paper) we take the cut-off scale 
to be A ~ Anv/ ^J\l — a^\. Since we are concerned with the phenomenology of single scalar 
production, the higher-order derivative operators are suppressed by powers of 0{m\/ Is?). As 
such, we are justified in neglecting such sub-leading effects in this paper. Non-derivative higher- 
dimensional operators will also exist in general. When the h field is not assumed to have a UV 
origin such as the SM Higgs, and is simply considered to be a singlet field [that need not neces- 
sarily transform under the nonlinearly realized SU(2)l x U(1) symmetry], the leading operators 
in the expansion in inverse powers of A appear at dimension five and are given by 

CU = -'-^hG^G^^'^ -'^hW^W'^^'^ -'-^hB.yB^'^ . (4) 

^ The symmetry assumptions we adopt by directly interpreting the data are minimal. In more involved scenarios, 
these assumptions can in principle be relaxed (see e.g. Ref. Il35l for a study with this aim that relaxes SU(2)c 
constraints). Note, however, that relaxing the assumption Cj y^j'^ = cy^j'^ significantly, can lead to conflict with 
precision constraints sensitive to SU(2)c and flavour violation. 
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Here gi, g2, gs are the weak hypercharge, SU(2)l and SU(3)c gauge couplings, respectively, and 
the different tensor fields are the corresponding field strengths with their associated Wilson coef- 
ficients Cj. The scale A corresponds to the mass scale of the lightest new state that is integrated 
out, which we assume is proximate to A. We will neglect operators originating from CP- violating 
sources due to the lack of any clear evidence of beyond the SM CP violation in lower energy pre- 
cision tests. Note that the operators in Cjjj^ can be further suppressed compared to the effects of 
a, c on (single) scalar production when the scalar field has specific UV origins. This is the case for 
example when his a PGB, see Ref. [|24]| for a detailed discussion. 

When one assumes that h is embedded into an SU(2)l doublet - if - as in the SM, the operators 
in first appear at dimension six, and the coefficients are suppressed by an extra factor ofv/A 
when considering single scalar production. In this case, the dimension six operator basis is also 
extended by the operator 

= -"^^j^HW^HB.^W^^r (5) 
For phenomenological purposes it is convenient to rotate to a basis for the operators given by 

Chd = -"-ifhG^^G^^^-'-^l^^hF.^F^^, (6) 

where F^^ is the electromagnetic field strength tensor and = cw + cb in the case of an SU(2)l 
singlet field, and = cw + cb — cwb if h is embedded into an SU(2)l doublet. 

In this manner, one can understand that the choices to retain the effects of higher dimensional 
operators (or not) in performing global fits introduces implicit UV dependence. In introducing 
higher dimensional operators, one is also explicitly assuming the existence of new states charged 
under at least a subgroup of the SM group. Alternatively, if NP is uncharged under the SM group 
but couples to the h'^ operator^, then it can impact /i-phenomenology by inducing an invisible h- 
width (when the scalar field takes a vacuum expectation value). This leads to the modification of 
the SM branching ratios for each decay into visible SM final states / via 

Bi{h ^ /) = (1 - Br,„,) X BisM{h -> /). (7) 

We will update our recent fit of Brj„„ to Higgs search data [1411 in a later section, and use this fit 
as a diagnostic tool to test aspects of our fit procedure. 



^ New physics of this form is sometimes referred to as coupling to the SM through the Higgs portal, see Refs. lf36H40l 
for some related discussion. 
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In general, the coefficients 'Q^inv ■ ■ ■ are arbitrary parameters subject to experimental 

constraints. The generic cases we consider are: 

• Composite/Pseudo-Goldstone Higgs/Dilaton scalar theories. 

In this case, a, c are free parameters in general, although they can be fixed in particular UV 
completions. If the scalar field is a Pseudo-Goldstone boson, it is also appropriate to neglect 
higher dimensional operators. We will fit to subsets of the parameters {a, c, BRinv, c^, Cg} 
in what follows. As this is a more general framework than the SM Higgs, we will use this 
EFT in assessing to what degree the SM Higgs hypothesis is consistent with the data or if 
deviations into this parameter space can give a substantially better fit. 

• The SM Higgs as a low-energy EFT. 

When the low-energy EFT is just the SM, the field h becomes part of a linear multiplet 



&3 = C2 = 0, and the only effect of NP is through non-renormalizable higher-order operators. 
The naturalness problem of the SM Higgs mass operator, and recent experimental hints of 
deviations in the observed properties of the (assumed) Higgs, motivates moderately heavy 
NP and the introduction of BSM parameters (c^, Cg, BRinv) as free parameters. We study 
the constraints on these parameters in detail in this paper. 

We will not attempt to relate the constraints obtained on the various parameters to any particular 
underlying model, other than the SM, in this paper. This choice is motivated by the current lack 
of other clear experimental evidence of BSM states to guide coherent model-building. The classes 
of models discussed above can be considered as motivating examples. 




(8) 



reducing Eq. ([2]) to the SM Higgs Lagrangian. In this case a = 6 = c = ^3 = ^4 = 1 and 



III. DATA TREATMENT 



A. Signal-Strength Data 



In this section we describe our method for globally fitting to the parameters discussed above, 
and incorporating the recently released 8 TeV data [[II, updated 7 TeV results from ATLAS [[42|. 
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the released 7 TeV CMS data P3l . and the recently reported Tevatron Higgs search results [|44l . 
This work builds on our previous fits flU |4T1. We only summarize the main details of the fit 
procedure and method here. Many subsidiary details of the fit procedure can be found in these 
reference works. 

We fit to the available Higgs signal strength data, 

_ Ei '^j^h X Br(/i i)]observed 
^ ^J^h X BT{h ^ Z)]SM ' 

for the production of a Higgs that decays into the visible channels i = 1- ■ ■ Nch, where Nch denotes 
the number of channels. The label j in the cross section, dj^h, is due to the fact that some final 
states are defined to only be summed over a subset of Higgs production processes j. The reported 
best fit value of a signal strength we denote by /i/. 
The global we construct is defined via 

x\f..) = y: ^^^^^ ■ (10) 
i=i ^» 

The covariance matrix has been taken to be diagonal with the square of the 1 a theory and ex- 
perimental errors added in quadrature for each observable, giving the error cTj in the equation 
above. Correlation coefficients are neglected as they are not supplied by the experimental col- 
laborations. For the experimental errors we use ± symmetric la errors on the reported fii. For 
theory predictions of the o-j^h and related errors, we use the numbers given on the webpage of 
the LHC Higgs Cross Section Working Group [|45l .^ The minimum (Xmin) determined, and 
the 68.2% (1 a), 95% (2 a), 99% (3 a) best fit regions are plotted as = xLn + with the 
appropriate cumulative distribution function (CDF) defining Ax^. 

We assume, as in Ref. [|2l[l3l, that the signal strength /ij in a given channel i follows a Gaussian 
distribution with the probability density function (pdf) given by 

pdf,(/i,,/i„cri)^e-(^-'^')'/(2.?)^ (11) 

with one-sigma error cTj, and best fit value /i,;. This is the case as long as the number of events is 
large, > (9(10) events, [[T3l . We normalize these pdf's to 1 in the interval (0, oo). 

In a simple counting experiment one has fii = {nobs,i — nbackg,i) , in terms of the observed numbers of events 
(nobs,i), the number of background events (nbackg,i) and the expected number of SM signal events (nf f^). 
^ These values have recently been updated for 7, 8 TeV and we use the updated numbers. Also note that BR(s s) is 
set to zero on this page but we use the latest version of HDECAY 14611 to add in BR(s s) to the quoted results. This 
has a negligible impact on the reported numbers through the modification of the total width. 
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In the framework of the SM, the predicted values of the /Xj are the same (equal to 1), and a 
universal signal strength modifier jj, can be defined and applied to all channels. By multiplying 
together the individual channel pdf 's (or the pdf 's of a single channel reported at two operating 
energies), we can also define a combined PDF for /i. This can be done for each separate experiment 
or for a global combination of all experiments. The combined PDF is also Gaussian and has 
combined fic and cxc values given approximately^ by 



We will use these relations to reconstruct the unreported 8TeV data from the reported 7 and 
7 + 8 TeV data.^ 

Armed with combined PDF(s), we can determine the 95% C.L. exclusion upper limits on the 
signal strength parameter /i (/x < HupL) OIEI. (We shall explain momentarily how to introduce 
an overall signal strength parameter in models with non-universal theory -predicted jii. In the dis- 
cussion that follows we are implicitly assuming that we are considering setting limits on combined 
signal strength parameters although we frame the discussion in terms of ji. A similar analysis can 
be carried through channel-by-channel instead of on the combined channels.) Such limits can be 
set on the combined signal strength parameter or for an individual channel's signal strength. 
With /i/s settling around unity with ai errors getting smaller and smaller due to increasing inte- 
grated luminosities, we can already start to set also lower limits on fi. The condition for such lower 
bounds, say at 95% C.L., to be meaningful is that the symmetric interval fi ± ^gs/t containing 95% 
of the integrated probability has /i — (Jgs/i > 0, in which case fJLdwL = A — ^gsA corresponds to the 
lower 95% C.L. bound on fj,. In this same case, we will take fi + ^gs/t > as the upper limit HupL- 
The conditions that define these limit are therefore 



^ This neglects (unsupplied) correlations and is therefore a rough approximation that should be taken with due cau- 
tion. An estimate of the accuracy of this procedure can be done by comparing quantities derived from such combi- 



is straightforward. These formulas are also easy to generalize in models with non-universal theory-predicted /i^. 
^ In version two of this paper we have incorporated a number of refinements in our treatment of the data, affecting 
the results that follow in the body of the paper. These refinements are now possible due to further information being 
released by the experimental collaborations after version one of this paper. See the Appendix for further details on 
the dataset now used. 
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nations vs. the experimental ones, which typically agree within 5-10 %. Generalizing ( [T2| to include correlations 
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and 



f^upL 



PDF(/i)rf;U 



erf 


f-^upL M 


1 + 


erf 


A 

V2a 



0.95/2 



(14) 



where erf(z) is the error function. However, when ft — 5<j^jl < 0, we shift the 95% C.L. interval to 
the asymmetric one (0, jjiupb) and revert to the upper limit definition [|2l[T3l^ 



^upL 



PDF(/i)d/i 



erf 


A 


— erf 


A- 


l-^iLpL 


1 + erf 


A 

y2a_ 





0.95 



(15) 



Next consider models that depart from the SM, like SM(a, c) with couplings of the Higgs to 
fermions and gauge bosons modified by the a, c factors as explained above. The predicted values 
of the yUj's will deviate from 1 in a channel-dependent way. In order to keep the same expected 
value of the signal- strengths for all channels, it is convenient to normalize the /i/s in Eq. ([9]) not 



to the SM signal expectation nf^^ but to the expectation in the model considered 



n 



SM(a,c) 



By 



doing this we can again use a universal signal-strength modifier /i, with expected value = 1, 
corresponding to the model being tested. On the other hand, with this change in the /i^ definition, 
the observed fii ± cTj values need to be rescaled by a factor n^f^ /n^^^"''^\ The individual pdf's are 
modified accordingly and can be combined in the usual way. 

To summarize, we can exclude (at 95 % C.L.) a given scenario not only if it predicts too many 
signal events which are not seen (/i > Hupl), but also if it predicts too few events, incompatible 
with the observed excesses associated with the reported discovery. The significance of such lower 
bounds will grow with more luminosity. Significances above 5a become possible (by definition) 
after discovery (which excludes the background hypothesis /i = 0). We will plot both of these 
bounds mapping the allowed fJ^upi, jJ^^wL ^^^^ '^he relevant parameter space through the dependence 
of Hi on the free parameters in our fit results. Figure |6| in Section |V} shows examples of such 
lower limits. 



The 95% C.L. interval extends down to when the ratio cr//t > 0.6, in which case there is no lower limit on and 



the upper limit is given by Eq. ( 15 i 
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B. Electro weak Precision Data 



We incorporate EWPD [|47l - |49ll by adding it directly to the measure in Eq. ( 10). When a is 



considered a free parameter, the shifts of the oblique parameters S and T are given by [|32ll 



6 TT 



AT 



3(1 - g^) 

8 TT C0S2 Qy^ V A 7 



(16) 



The numerical coefficient is determined from the logarithmic large-m/j dependence of S, T given 
in Ref. [|48i.^ As for EWPD, recent updates to the measurement of m^y at the Tevatron [[SOllSTIl 
have refined the world average [|52l . and have significantly reduced the quoted error. Incorporat- 
ing"' these new measurements we use [|53l 



5 = 0.00 ±0.10, 



T = 0.02 ±0.11, 



while the matrix of correlation coefficients is given by 



C 



1 0.89 
0.89 1 
-0.55 -0.80 



-0.55 
-0.80 
1 



V = 0.03 ±0.09 



(17) 



(18) 



Here we have assumed = 125 GeV, as corrections for shifting these results by a few GeV 
are negligible. There is a strong preference for a ~ 1 in the global fit when EWPD is used, and 
the constraints on the scalar field can be directly associated with EWPD bounds. Note that the 
slight preference for a > 1 in the best fit region when EWPD is taken into account is subject to 
uncertainties in cut-off scale effects. Although the shift in the best fit point is of interest as a probe 
of possible new physics, we cannot clearly disentangle such a hint from cut-off scale effects. 



IV. FIT RESULTS 



A. Status of the Higgs hypothesis 

The excess of events of ^ 5a significance reported by ATLAS and CMS peaks, as a function 
of the scalar mass, at slightly different values: 126.5 GeV for ATLAS and 125 GeV for CMS. 

^ Here we have introduced an Euclidean momentum cut-off scale A. The degree to which A properly captures the 
UV regularization of S and T is model-dependent. We assume that directly treating this cut-off scale as a proxy 
for a heavy mass scale integrated out is a good approximation, i.e. that further arbitrary parameters rescaling the 
cut-off scale terms need not be introduced. 
We thank J. Erler for kindly providing these EWPD results. 
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This difference can be attributed at this stage of the search to statistical fluctuations in the data. 
Monte-Carlo studies ll54ll indicate that such effects can shift the observed maximum signal strength 
compared to the true signal strength maximum by ~ 2 — 3 GeV. Due to this, a fit that combines 
data from different experiments at the same mass value might be biased and not necessarily better 
than using data at slightly different masses. As much more detailed data is available to us at the 
mass peaks, in our global fits we will use all available fXi taken at ruh = 125 GeV for the CMS and 
the Tevatron (which has a ~ 3cr excess over a wider region of masses), and at rrih = 126.5 GeV 
for ATLAS (see Fig.|7|in the Appendix)." 

In order to assess the degree of consistency of the data with the SM hypothesis we first con- 
sider the effective Lagrangian given by Cejf and assume that higher dimensional operators are 
sufficiently suppressed so that C/yj Cg C3.n be neglected. We then perform a two-parameter x^(a, c) 
fit and examine the A for the SM point (a, c) = (1,1) compared with the best fit point. This 
defines a C.L. corresponding to the deviation of the SM hypothesis compared to the best fit point. 
The result is shown in Figure [T] (left). We also show in Fig. [T] (right) the best fit regions when 
EWPD is added to the global measure. Notice the dramatic reduction in the size of the best fit 
region along the a-parameter, which is forced to lie close to 1. These results visually summarize 
the current experimental status of establishing the Higgs hypothesis. 

When EWPD is not used, the SM Higgs hypothesis of (a, c) = (1, 1) is < 2 (C.L. of 0.88) 
away from the best fit point, which sits at (a, c) = (1.1, 0.68). Note that here and in the following 
discussion we are choosing to round the C.L. and the best fit points. This is due to the preliminary 
nature of the 7 + 8 TeV data and is not limited directly to this accuracy due to the fit procedure 
we have adopted. The C.L. of the SM hypothesis in the combined data is consistent with our past 
results at 7 TeV With the updated data released and incorporated in our fit since version one of 
this paper the parameter space that has the global minimum has changed from c < (with initial 
ICHEP data) to c > (with post-ICHEP updates).'^ The existence of the non negligible parameter 
space with c < is easy to understand. Due to the interference term in the /i — t- 7 7 decay width 
which is oc —ac, a negative c allows a relatively larger excess in 77 events due to constructive 

" Recent updates from the experimental collaborations that have appeared since version one of this paper have re- 
leased more information at ruh = 125.5 GeV for CMS [55 1 and at ruh = 126 GeV for ATLAS |l56l. However, 
as the amount of information relevant to the fits we will perform released to date is still greater at the mass scales 
rrih = 125 GeV for CMS and nih — 126.5 GeV for ATLAS we retain these mass choices in our fit. 
The lack of data released for these subcategories to date at rrih = 126 GeV is the primary reason we retain the use 
of = 126.5 GeV for ATLAS data. 



12 



7&8 TeV LHC data & Tevatron 7&8 TeV LHC data & Tevatron + EWPD 




0.6 0.8 1.0 1.2 1.4 0.95 1.00 1.05 1.10 

a a 



FIG. 1: Global fit results in the (a, c) plane for all reported best fit values given by ATLAS and CMS, left 
(right) without EWPD (with EWPD). In both plots we take ruh = 125 GeV for the Tevatron and CMS7/8 
and rrih = 126.5 GeV for ATLAS7/8. The green, yellow, gray regions corresponds to the allowed 1, 2, 3 cr 
spaces for a two parameter fit. The best fit point in each region is also labeled with a point. The thicker 
point indicates the one with the smaller Xmin- 

interference between the top and W boson loops. When EWPD is used as in Figure [T] (right) we 
find that the SM is similarly residing at ~ 2 cr (C.L. of 0.93) away from the best fit point which is 
now (a, c) = (1.0, 0.67) and the best fit region where c > now has a (significantly) lower global 
minimum. The minima are no longer as degenerate with the addition of the most recent ATLAS 
data, Ax^(mmi,mm2) ~ 4. 

In view of the different masses of the signal- strength peaks in the various experiments (which 
can be due to the statistical effects mentioned above) and of the subtleties we have neglected in 
properly combining the results of these different experiments, it is also of interest to perform the 
fit in the (a,c) space for each experiment individually. We show these results in Figure [2] The 
CMS experiment has the SM point residing about ~ 2a from the best fit point, with the C.L. of 
the SM case compared to the best fit point at 93%. For ATLAS, the SM point is now at a C.L. of 
41%, within the ~ la region. The Tevatron results have the SM point within the la region with a 
C.L. of the SM case (compared to the best fit point) of 50%. 

The allowed fit region for CMS can be compared to the recently presented public results [HI, 
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7&8 TeV LHC data CMS (mh=125) 7&8 TeV LHC data ATLAS(mli=126.5) Tevatron Data, = 125 GeV 
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FIG. 2: Best fit regions (at 68%, 95% and 99.9% C.L.) in the (a, c) plane for a fit to all reported 
signal-strength values given by ATLAS (m/j = 126.5 GeV), CMS (nih = 125 GeV) and the Tevatron 
(mfi = 125 GeV) collaborations individually. We plot the same best fit contours over the same domain of 
parameter space to allow a direct comparison amongst experimental results. The significant change in the 
ATLAS results from version one of this paper is due to the use of the ATLAS diphoton data broken into 
subcategories. 

which restrict the fit to the region c > (physically different from the region c < 0) As the 
absolute minimum of the lies in fact in the discarded region, the shape and size of the 68% and 
95% C.L. regions presented by CMS differ from the ones that we obtain in figure |2] If we also 
restrict our parameter space to positive c we have checked that we get excellent agreement with 
the CMS result. Note however that there is no valid reason to discard a priori the negative c region 
which offers in fact the interesting possibility of giving a good fit to the data by an enhancement 
of the h'j'j coupling through constructive interference of the top and W loops. 

B. BSM Implications 

1. Implications for an invisible width 

The results of the last section can be interpreted as (partial) evidence in support of the SM 
Higgs hypothesis. Assuming then that the observed boson is the SM Higgs, we can study possible 



We can always set a > by a redefinition of the h field. The sign of c could in principle be changed also by rotating 
the fermion fields, but this would affect in the same way the fermion mass so that the sign of c/m/ is in fact fixed 
and physically meaningful. 
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FIG. 3: Global fit to Bvinvfor the SMHiggs using only CMS data (left figure) for nih = 125 GeV with two 
methods as a check of our fit procedure. Blue curve - global signal strength based fit, Red curve - individual 
channel fit. See text for further explanation. The middle figure shows the distribution developed from the 
combined best fit fic supplied by the four experiments, including the 7 and 8 TeV LHC results for two mass 
values. In the right figure we show the discovery potential for Brj„^, updating a result from the analysis in 
Ref. / |47I/ with the new global signal strength data. 

deviations of its properties due to BSM effects. The simplest extension of the SM in terms of 
new parameters is perhaps the case where one only introduces a Higgs invisible branching ratio, 
Bimt,. Such an extension is common in many BSM scenarios, e.g. when new physics couples 
through the Higgs portal, and new states exist that are uncharged under the SM group. Fitting to 
Bijnt, allows one to fit to the global combined signal strength values supplied by the experimental 
collaborations. The values we use to perform the fit are given in Table |lj 



Experiments 


fic, ruh = 125 


fJc, nih = 125 


flc, nih = 126.5 


cjc, mh = 126.5 


CMS [7&8 TeV] HJ 


0.80 


0.20 


0.67 


0.19 


ATLAS [7&8 TeV] ffl 


1.12 


0.27 


1.24 


0.26 


ATLAS [7&8 TeV] (& iiww) ESl 


1.32 


0.29 


1.37 


0.27 


CDF&D0 1441 


1.35 


0.59 


1.38 


0.60 



TABLE L Combined signal strengths fic and errors cjcfrom ATLAS, CMS and the Tevatron collaborations. 
Here we quote it symmetric la errors. 



The left plot of FigjS] shows the result of extracting Brj^^, using our fit approach using two 
different methods. This allows an important cross check of our procedure and results. We have 
employed in our fit two approximations that require further justification. Using the assumption of 
Gaussian PDF's, we have extracted the 8 TeV data from the known 7 TeV data and the released 
7 + 8 TeV data. In doing so we have neglected (unsupphed) correlations. Further, in performing 
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our fits we have neglected correlation coefficients for the /ij. The left plot of FigjS] shows in blue 
the result of fitting directly to the combined CMS data - which do take into account correlations. In 
red we show, for comparison, the fit results when we use our procedure to extract the 8 TeV data, 
and fit to the individual /x/s. The best fit point and the 95% C.L. regions are in good agreement, 
supporting the estimated accuracy of our results of ~ 5 — 10%. 

The middle plot of Fig js] shows instead the resulting distributions for Brj„^ for two different 
mass values = 125, 126.5 GeV when we combine the results from the three experiments and 
fit to the supplied signal strength parameters flc, o^c- This is the most accurate analysis on Brj„^ 
that we can perform with the released data (more accurate than a fit on the individual channels). 
This is primarily due to the experimental correlation effects that are incorporated in the fic, and 
larger number of channels that are incorporated in the experimental likelihoods used to construct 
the combined signal strength parameter used. We find the 95% C.L. regions Bimv < 0.37(0.40) 
for nih = 125(126.5) GeV. These limits can be used to constrain many models that predict an 
invisible Higgs width. For recent related work, see Ref. lfT5ll4Tll58l . 

The right plot of Fig j3] (which updates a similar analysis in Ref. 011) shows the current status of 
the quantity (1 — fic) (with one sigma error band), interpreted as an invisible Higgs branching ratio, 
as a function of the Higgs mass. This result is also based on the combined signal strengths supplied 
by the experimental collaborations. With the current errors, no significant statement can be made 
about the possible nonzero central value of Br iny. The shaded areas indicate the 2a range of 
possible fluctuations (of the background in the upper region at large Bij^^; of the SM Higgs signal 
in the lower region at small Bij^j,) that could be miss-interpreted as an invisible Higgs width. This 
plot shows that the current nonzero central values of Brj„^, at nih ~ 125 GeV (consistent with the 
fit results) are perfectly compatible with downward fluctuations of the SM Higgs signal. 

The main difference between the results presented in Figj3] with respect to the first version of 
this paper comes from the use of the ATLAS combined signal strength reported in Ref. [15611 . that 
incorporates now the large 8 Tev ATLAS WW 153 signal strength fiwwimh = 125) = 2.1+°;^. 
One finds the 95% C.L. limits Br^^^ < 0.32 (0.37) for ruh = 125(126.5) GeV when this signal 
strength is added to the global data set. 
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2. Implications for Cg,c^ 

One can also infer the current experimental bounds on the BSM parameters Cg, c^. We expect 
that these operators arise at the loop level, so we rescale the Wilson coefficients as Cj = Cj/ (IGtt^) 
for j = (7,7. Using the results of Ref. [|59ll . the effects of these operators are incorporated as 
rescaling factors used in the fit and given by 



Ra = 



g - SM 
^gg^h 



1 V^Cg 



7 — p5M 

/i— >7 7 



0.75ci - (0.05 + 0.072)ci, A2 

1+ ^ 



(19) 



(20) 



4(2.07a - 0.44q + (0.01 + 0.01 i)cb) 
Here we have used rrit = 172.5 GeV, = 4.75 GeV, nih = 125 GeV and a,(172.5 GeV) = 
0.107995. When the Higgs has SM renormalizable couplings, then a = q = Cf, = 1. We have 
retained the two-loop QCD correction to the SM matching of the hG^j^G^^'^ operator in the 
nit ^ 00 limit in these numerical coefficients. The operators in C%j^ also affect Br(/;, — )• 7Z), 
but the effects in our fit can be neglected. See Ref. [[411 for further details and discussion on our 
fitting procedure to these higher dimensional operators. 

We show in Figure|4]the results of fitting to c^, using the current dataset when one assumes 
the SM values a = c = 1 (left); when one fits to and subsequently marginalizes over a, c 

(middle); and finally, when one fits to Cg,c^, BRj^t, and marginalizes over BRm^,. These results 
summarize the preference in the current dataset for including higher dimensional operators when 
the scalar is not assumed to be the Higgs. 

We also show in Fig. [5] (left) the fit to c^, c^, a, c, marginalizing over the higher dimensional 
operators. This case shows the preference for a, c even when more massive states are integrated out 
of a composite scalar theory (for example) contributing to Cg, c^. We see that the SM hypothesis is 
significantly improved in its consistency with the dataset in the context of NP of this form. This is 
particularly relevant as it directly shows what the data tells about a, the key parameter that probes 
the involvement of h in electroweak symmetry breaking. A final plot in Fig. |5] (right) shows the 
allowed space when q is varied independently from the remaining fermion couplings which are 
treated with a universal rescaling. 

The fit to higher dimensional operators shows a preference for a BSM contribution to c^, which 
deserves some comment. First, one can study the model independence of this preference for 
an enhancement of by constructing a one dimensional distribution, marginalizing over the 
unknown operator Cg when the SM values a = c = 1 are assumed, or over (c^, a, c) when one is 
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7&8 TeV LHC data & Tevatron 



(a,c) marginalized 



7&8 TeV LHC data & Tevatron BRi^, marginalized 



cJA- 





-30 -20 -10 10 20 



FIG. 4: (Left) Results of fitting to Cg,c^ wfien tfie SM is assumed. Note that in these results we have 
assumed that Cg, are real. (Middle) Results of fitting to Cp, c-^, a, c and marginalizing over (a,c) subject 
to the constraint a > and < c < 3. (Right) Results of fitting to Cg, c^, BRj„„ and marginalizing over 



V Cg/A , V Cy/A marginalized 



7&8 TeV LHC data & Tevatron 



c 







4 





2.5 F 




FIG. 5: (Left) Results of fitting to Cg,c^,a,c and marginalizing over the higher dimensional operators. 
(Right) Allowed space when a is varied independently. See text for more details. 

considering non-SM scalar scenarios. Performing this exercise we find that a preference for c-y ^ 
only exists in certain cases where implicit UV assumptions are adopted due to the parameters used 
to fit the data. One can make a number of observations regarding enhanced /i — )► 77 event rates. 
The coupling to F^i^ can come about due to cwb- If this is the case, one can bound the allowed 
enhancement of /i^^ due to related EWPD constraints. Using the relation [|59] - l6n 

^,2 



A2 



Cwb = -2n S, 



(21) 
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one finds the current experimental constraint \v'^ cwb/^^\ ~ for an Na deviation consid- 

ered in the EWPD parameter S, assuming rrih — 125 GeV. Using //^^ Rg Rj, this leads to 
the bound < 1.2 (1.4) for a one (two) sigma deviation in the S parameter when Cg = 0. Re- 
stricting the unknown higher dimensional operators to have a global measure in the 2a allowed 
region (Ax^(cg, Cj) < 6.18), one can maximize the enhancement of fx^y^ considering the related 
constraints on c^^b. One finds that an enhancement of 1.4 is possible when a one sigma deviation 
in the S parameter is also allowed. The Cj Wilson coefficient could also come about due to cb,cw 
or a combination of the two. In this case, the EWPD constraint does not directly apply, and an 
enhancement of //^^ by a factor of < 1.7 is still allowed when maximizing the contribution subject 
to the constraint Ax^{cg, c^,) < 2.3. 

This analysis has assumed no relationship between the Wilson coefficients c-y, Cg, which is fixed 
in any particular UV completion. It is interesting to note that in general when integrating out a 
single BSM field one expects a strong relationship between the Wilson coefficients, with identical 
loop functions in many cases, the only differences in the matching onto the Wilson coefficients is 
due to the relative charges of the new states under the SM subgroups SU(3)c and U(l)em- 

V. STUDIES OF CONSISTENCY AND TENSION WITHIN THE DATASET 

The consistency between the search results (fii ± cTj) from different channels can be quantified 
as follows. For each channel i we construct its Gaussian approximation to the pdf of the signal 
strength, pdf j(/i) = pdf (//, /tj, ctj), which we can contrast with the full combined PDF (either in a 
given experiment or for the overall combination of all data), PDF(//) = /tc, Cc). For each 

single channel we calculate its p-value with respect to the global combined PDF, i.e., pic is the 
p- value for fi-i assuming the full PDF(/x). We can also define the p- value of the global jic with 
respect to the individual pdfj(yu), which we will denote by p^. In order to treat properly the p- 
values for cases with fii < 0, we will normalize the pdf 's in this section so that they give 1 when 
integrated over the whole interval /i e (—00,00). We then have, for fii < jlc. 




(22) 
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7&8 TeV CMS data (mh= 125 GeV) 7&8 TeV ATLAS data (mh= 126.5 GeV) 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



a a 



Tevatron data (mh= 125 GeV) 7&8 TeV LHC & Tevatron data 




0.0 0.5 1.0 1.5 2.0 2.5 0.6 0.8 1.0 1.2 1.4 



a a 

FIG. 6: Exclusion limits on the SM(a, c) parameter space derived from a) 95% C.L. upper limits on the 
signal strength parameter (i < fiupL (gray shaded regions beyond the black lines towards the right), b) 
95% C.L. lower limits fi > ^dwL (gray shaded regions enclosed by the black lines close to the origin) and 
c) "tension limits" (purple dotted regions delimited by purple straight lines) from the presence of search 
channels inconsistent with the rest of the dataset at 2 a. The most powerful channels in setting such hmits 
are listed in the text. For comparison, the best-fit regions discussed in previous sections are also shown. The 
fermiophobic scenario along c = is excluded. 
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and, for /tj > /ic, 



Pic = I PDF(/i)(i/i = -erfc 



2 V y2^T. 



Pci = / pdfj(/i)d/i = ^erfc ( ) , (23) 



— oo 



2 V V2f^ 

where erfc is the complementary error function, erfc(z) = 1 — erf(z). We will say channel i is 
in tension with the rest of the data if pic and pd are both very small. For a given critical p-value 
Pn, corresponding to N standard deviations, channel i is not consistent with the combined dataset 
at {1 — pn) C.L. if Pic, Pci < Pn and the model can be excluded based on that disagreement. The 
consistency condition reduces simply to 

\fii - fic\<N max[(Ti, (Tc] . (24) 

We will choose = 2 in our discussion. 

When this test for consistency is applied to the {fii ± cij) dataset, interpreted as coming from 
a SM Higgs signal, we find tension at this 2o level for two ATLAS 77 channels at 7 TeV, those 
labelled URhPTh and CChPTt (the two outliers easily identified in Fig. |7]). In extensions of the 
SM, like the two-parameter scenario SM(a, c) we have discussed in previous sections, the rescaling 



of the different channels (which affects jli ± (Xj as explained in section III A I can introduce very 
significant distortions in the pdfs and cause too large tensions for some other channels. Such 
regions of parameter space can therefore be excluded on this basis. Figure [6] illustrates this by 
showing, on addition to the best fit regions and the 95% C.L. exclusion regions derived from upper 
and lower limits on the signal strength, the regions in (a, c) space that would be excluded due to 
more than 2-cr tension in some search channel (purple shaded regions delimited by straight lines). 
We show such limits both experiment by experiment and for the combined result. Typically, in 
the excluded regions several channels at the same time cause the exclusion. We show in each case 
all the region that is excluded by at least one channel (with the exception of the ATLAS case, 
see below). In each case, the channels that have bigger exclusion power are the following. For 
Tevatron, the hh channel; for CMS, at 7 TeV, 77^^ and at 8 TeV, 773 and rr; for ATLAS, the 
two 77 channels mentioned above, URhPTh and CChPTt at 7 TeV, cause 2a tension in all the 
region of parameter space shown and we do not mark this area. Besides these channels, there is 
2 a tension also from 'yycciPTt at 7 TeV and from WW -> lluu at 8 TeV. The fact that the tension 
limits are straight lines passing through the origin is due to the fact that any common rescaling 
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of a and c (that leaves c/a invariant) also leaves invariant the tension associated to any channel 
as the latter are functions of ratios (/ij — fJ,c)/o'i,c, and such ratios are also functions of c/a. We 
see from Figure [6[ that fermiophobic scenarios, along the axis c = 0, are excluded at this level of 
confidence. 

This channel-tension analysis is clearly related to the study we have also performed and 
the tension exclusion limits tend to exclude regions of parameter space that give a bad fit to the 
data. However, this approach seems to be more powerful in being able to exclude definitely some 
models. We also see that there are regions of parameter space that are not excluded by the con- 
ventional upper (or lower) limits imposed on the signal strength parameter, yet can be excluded by 
the tension exclusion limit. This is simply due to the fact that the dataset can contain two widely 
separate channels and still give a combined PDF that respects the upper limit, which is not able 
to probe in such cases the internal inconsistency of the individual channels. We conclude that this 
type of analysis offers a complementary tool in testing model performance. 

VI. CONCLUSIONS 

We have studied the evidence of a scalar field that has been discovered by the LHC collabo- 
rations using an effective field theory framework. We have discerned what can be inferred about 
the properties of the scalar field at this time using joint fits to available datasets. We have 
also developed and applied new techniques that allow one to exclude model classes that are in 
tension with the data through violating upper or lower C.L. bounds, or through introducing exces- 
sive tension into the signal strength parameter datasets. At this time, according to our fit method, 
and using publicly available data, the SM Higgs hypothesis is consistent with the global dataset 
compared to the best fit point at the level of < 2cr. 

Appendix A: Data Used 

Our approach to the presented data is as follows. As the relative weights of the various contri- 
butions to the inclusive Higgs production processes depend on the operating energy of the LHC, 
we need separate information on 7 and 8 TeV data to perform our fit as the relative fraction of 
events at each operating energy depends on the unknown parameters a,c,Cg,Cy}^ Most of the 

^'^ Fits to BRinv do not suffer this problem as this parameter affects universally all channels. 
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public results (T\ are now presented separately in 7, 8 TeV signal strengths. When this information 
is available we use directly this data. 

For those cases in which only the combined 7 + 8 TeV signal strengths are available (currently 



CMS combined firr, P'WW used in the fit) in addition to the 7 TeV results, we make use of Eq. (12) 
to reconstruct the unreported 8 TeV data. Note that this relies in the use of Gaussian approxima- 
tions to the PDF's to describe the data (and signal strengths) and should be increasingly accurate 
as the total number of events increases. This can be done without knowing directly the experimen- 
tal likelihood function in the limit that correlations are neglected, which is an assumption we are 
already forced to adopt as this information is not supplied by the experimental collaborations. We 
show in Fig. |7] the resulting reconstructed data. 

An interesting check of our approach is to use a subset of the provided subclass signal strengths 
to reproduce a reported combined signal strength. This exercise can be carried out, for example, 
on the supplied vector boson fusion tagging b b signal strength and tth signal strength (that uses 
h ^ bb) and comparing to the presented combined h ^ bb signal strength. We have carried 
out this procedure for the presented CMS data and find good agreement with the reported results, 
within our estimate of 5 — 10% error introduced due to a lack of correlations. 

We have updated our fit from vl of this paper to include the following information that has 
been released since the first version appeared on the archive. The full subclasses of 7 7 events are 
now available from ATLAS and CMS at 7, 8 TeV. Also, the production channel composition of 
the 77 subclasses have been supplied [|63ll64ll . We have modified our fit to use this information 
consistently with our rescaling procedure. This procedure replaces our utilization of estimated gg 
contamination in the pp — > 7 7 j j signal events. Also note that despite the cuts of the "tight" and 
"loose" dijet channels of CMS indicating they are not mutually exclusive event classes, the CMS 
collaboration vetoes an event appearing in the "loose" pp — )■ 77 jj sample if it passes the "tight" 
cuts [65|. As such, these event classes can both be included in the global fit we perform.'^ Fi- 
nally we note that the 7, 8 TeV signal strengths for the b b and ZZ signal strengths have now been 
supphed and are incorporated in our fit. Comparing our extracted ZZ 8 TeV result to the experi- 



We however still thank V. Sanz for kindly providing the contamination coefficients for vl appropriate for our 
past procedure.The contamination in the first version of this paper was reported with a typographical error, the 
correct contaminations are e = 0.032 for the 7 TeV dijet tagged diphoton signal, e ~ 0.023 for the 8 TeV "tight" 
dijet tagged diphoton data and e = 0.039 for the 8 TeV "loose" dijet tagged diphoton data 1621 . Here e is the 
contamination of the — > 7 7 jj signals due to gg Higgs production events, when e is defined such that the rate is 
given by {ecTgg^h + cTjjh) x Br{h 77). 
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mentally supplied number we find agreement within the estimated accuracy of our procedure. For 
the b b CMS signal strength we note that the 7 TeV signal strength recently reported in Ref. Il66ll 
differs from the previously public 7 TeV signal strength. 



All Data (48 channels) 



WW->llvv ; 0.32±0.83 
bb : l.97±0.71 
yy : 3.62±2.76 



TT : 0.45±1.8 
bb : 0.46+2.18 
ZZ->41 : 1.11 ±0.92 
WW->llvv : 0.45±0.62 
TTuciPTi ■''■54±l-46 

r>'DRiPT,-2-52±1.69 
rXmiFT. ■■ 10-42±3.68 
rrcam ■ 6 08±2.63 
>Tca,Pr, : -4.36±1.78 
rJ-cR,™ : 2.74±2.a3 
rr™„, : -1.63±2.87 
yy^^ : 0.35±3.61 
ry,, :2.71±1.9 



ATLAS 7 TeV 



ZZ^41 ; 1.03±0.64 
WW->llvv : 1.85±0.65 
ryuciPT, : l-04±l-3 

7yj^p„:0.95±1.39 

TYcavr, ■ 3.48±1.95 
rrccKT, : 3.57±2.59 
rW,:0-47±1.74 
yrcRhFr,:0-4±2.1 
TTct : 5.55±3.3 
yy.. : 2.58±1.73 



ATLAS 8 TeV 



TT : 0.62±1.13 
bb : 0.59±1.17 
ZZ^41 : 0.64±0.57 
WW->nw : 0.38±0.56 
•)T„:3.16±1.81 
77, : 0.66±0.95 
77, : 0.74±1.I6 
77, : l.53±1.6 
77,, : 4.2 1 ±2.05 
TT : -0.72±0.97 
bb : 0.41 ±0.94 
ZZ->41 : 0.79±0.56 
WW->llvv : 0.98±0.71 
r7„ : 1.45±1.25 
77, : 1.5±1.04 
772 • ('■95±1.16 
77, : 3.76±1.77 



jj tigiii 



1.32±1.58 
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FIG. 7: Pictorial presentation of the data used in the fits to sub-channels. Blue: reported data at 7 TeV. 
Red: reported 8 TeV data, or reconstructed 8 TeV data. 
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